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Abstract Mirror isobars 7 Li and 7 Be are investigated in a dicluster model. The magnetic dipole moments 
and the magnetic dipole response to the continuum are calculated in this framework. The magnetic contri- 
bution is found to be small with respect to electric dipole and quadrupole excitations even at astrophysical 
energies, at a variance with the case of deuteron. Energy weighted molecular sum rules are evaluated 
and a formula for the molecular magnetic dipole sum rule is found which matches the numerical calcula- 
tions. Cross-sections for photo-dissociation and radiative capture as well as the S-factor for reactions of 
astrophysical significance are calculated with good agreement with known experimental data. 
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1 Dicluster description of 7 Be and 7 Li 



It is well known that the properties of A = 7 mass nu- 
clei may be effectively described in terms of a dicluster 
model: two (inert Iclusters in interaction with each other 
CD; 0; i; II S SB El)- A long-standing description of 
these nuclei with such a dicluster picture has achieved ex- 
cellent results and a continued interest because it catches 
the essential physics within a very simple, yet powerful, 
model. Such a simple nuclear structure scheme may then 
used in reaction calculations (0). The need to go beyond 
this level of description may be justified only with the aim 
of taking into account more refined considerations on the 
role of the Pauli principle or on the occurrence of com- 
plex phenomena like polarization of the clusters, and not 
because the model itself does not reproduce static and dy- 
namic properties of A = 7 nuclei to a good degree of accu- 
racy. For example ab initio shell model calculations using 
a computationally-heavy Monte Carlo method have suc- 
cessfully reproduced many properties of these and other 
light nuclei starting from a nucleon-nuclcon interaction 
(flOl ). By computing ground state overlap functions, these 
calculations have shown that 7 Be has a spectroscopic fac- 
tor of almost 1 in the a+ 3 He configuration and, similarly, 
7 Li has a large spectroscopic factor in the a + t channel. 
A number of other approaches, that we will not discuss 
here, such as Resonating Group Methods or Generator 
Coordinates Methods, are available to various degrees of 
complication and success. In addition we must mention 
a very recent study by Canton and Levchuk on the low- 
energy capture reactions (fllh using the Multichannel Al- 
gebraic Scattering approach that has been applied to the 



3 He(ci!, 7) 7 Be reaction at astrophysical energies achieving 
a good agreement with S-factors data. 

We have shown in recent papers (d) that the dicluster 
model may be profitably used to calculate a number of 
static and dynamic properties of 7 Li, especially in con- 
nection with excitations to the continuum in break-up 
processes. We apply here the same model to 7 Be and in 
addition we investigate the magnetic excitations to the 
continuum that were neglected in previous works. It is in 
fact important to give proper estimates of their contribu- 
tions to reaction cross-sections, especially at energies of 
astrophysical interest. These two isospin mirrors not only 
display very similar ground state properties, but also man- 
ifest a similar behaviour as long as their response to the 
continuum is considered. Nowadays it is becoming more 
and more evident that to properly describe structure and 
reactions of loosely bound systems, e.g. nuclei along the 
drip lines, it is essential to take into account the coupling 
to the continuum and the break up channel. This is rela- 
tively simple in the case of one-particle halo nuclei. In the 
case of more complex weakly-bound nuclei such as 7 Li or 
7 Be, the cluster picture offers the possibility of describing 
the ground and a few excited states in terms formally anal- 
ogous to that of a single-particle picture, with the caveat 
that the single particle is not just a neutron or a proton, 
but a composite one, like for instance a triton. 

Of course, as in any theory driven by phenomenological 
considerations, a certain degree of approximation comes 
into play, and one needs to consider the pros and cons of 
the approach. We mention therefore that considering the 
clusters as spherical and elementary, but not point-like, 
has some implications: i) one neglects possible polarization 
effect, which may come into play because of the proximity 
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of the two clusters and ii) one neglects nucleon-exchange 
effects between the two clusters. We shall prove that these 
limitations are not severe (see also (Q)) and the model 
works rather well for the A = 7 nuclei. 

Our intercluster potential is assumed to be the sum 
of a nuclear Woods-Saxon potential, a Coulomb potential 
and a spin-orbit interaction: 



results for the photo-dissociation and radiative capture 
reactions. In Section IV we present the S-factors for the 
two reactions 3 He(a, 7) 7 Be and 3 H(a, 7) 7 Li in comparison 
with known experimental measurements. The conclusions 
are followed by an appendix where a few useful formulas 
for dicluster nuclei are derived. 



V Ai -a 2 = V cou i{r) + V W s(r) + Vi. s (r) 



(1 . 1) 



with V\. s (r) = Vi. s l ■ s^- dVw d ^ r ^ . Here r is the same ra- 
dius parameter that enters into the parametrization of the 
Woods-Saxon potential, t and s are the orbital angular 
momentum and spin of the 3 H( 3 He) cluster in 7 Li( 7 Be) 
respectively. The Coulomb potential has the standard 1 Jr 
behaviour at large distances (point charge), but at small 
distances we use the Coulomb potential generated by a 
uniform spherical charge distribution. 

Note that it is not possible to identify one of the two 
clusters as a core and the other as a valence particle, since 
the two fragments have comparable masses and charges. 
In the spin-orbit interaction, r is the radius of the Woods- 
Saxon potential, 1 refers to the relative angular momen- 
tum, and s refers to the spin of the triton (helium-3) which 
can be seen as a "heavy" single-particle for 7 Li ( 7 Be). 

The A = 7 isobars are particularly challenging. They 
are weakly bound, but nevertheless they have one excited 
bound state. The low-lying response presents (as we will 
show) non-resonant as well as resonant peaks with compa- 
rable intensities. Aside from the importance in pure nu- 
clear physics, the break-up, photo-dissociation and cap- 
ture processes for these nuclei play an important role for 
astrophysics and for applicative purposes (e.g. lithium is 
used as coolant in nuclear reactors). We apply, in par- 
ticular, our dicluster model to the calculation of astro- 
physical S-factors for the reaction 3 He(a,7) 7 Be (and also 
3 H(a,7) 7 Li). Recently the LUNA collaboration ([H) has 
undergone a successful experimental campaign aimed at 
measuring the cross-section for this reaction which is thought 
to i) strongly influence the hydrogen-burning process in 
the Sun and consequently the solar neutrino production 
and ii) determine the primordial abundance of light ele- 
ments which influences the determination of the baryon- 
to-photon ratio of the Universe from comparison with nu- 
cleosynthesis. In both these fields the determination of 
this quantity is crucial and this fact justifies even more 
the need for simple, but yet powerful, theories that could 
explain and reproduce the data. For a review of experi- 
ments in this direction see Ref. (QJ). 

The paper is organized as follows: the present section 
contains a summary of the method and deals with the 
parameters of the cluster-cluster potential and with the 
electromagnetic properties of the 7 Li and 7 Be systems. 
Static properties and transition rates among bound states 
are compared with data whenever possible. Section II con- 
tains a digression about molecular sum rules, their evalua- 
tion in the present case and the derivation of a new energy 
weighted sum rule for magnetic dipole transitions, inspired 
to the Kurath single-particle sum rule ((2ll :l22l). Section III 
deals with the application to astrophysics and presents our 



1.1 7 Be 

We consider 7 Be as formed by a and 3 He. The parameters 
of the cluster-cluster Woods-Saxon potential are Vws = 
-73.851 MeV for the depth, 1.60 fm and 0.48 fm for its 
radius and diffusivity. The magnitude of the spin-orbit po- 
tential, Ui. s = 1.275 MeV, is chosen in order to reproduce 
the splitting of the two bound eigenstates, respectively at 
— 1.58 and —0.98 MeV. The shape of the potential has 
been tuned in order to obtain an overall agreement with 
various observables, like the charge and mass radii. The 
a-particle has an intrinsic spin equal to zero, 3 He has spin 
1/2 and the intercluster orbital angular momentum for 
the ground-state is I = 1 : this can either be postulated on 
the basis of the Wildermuth connection (2) or it can be 
inferred phenomenologically, for example, from the mag- 
netic moment for a dicluster nucleus. Using the intrinsic 
magnetic moment, masses and charges of the two clusters 
in formula (|5 . 20| for the total magnetic moment, one 
can see from Table ([I} that the value that gives the best 
agreement is clearly 1=1 (which also gives the correct 
parity) . Due to spin-orbit interaction the ground state has 
j = 3/2 _ , while the first excited state has j = l/2~ . 





s 


P 


d 


exp. 


'Li 

7 Be 


2.97 
-2.12 


3.37 
-1.53 


3.78 
-0.93 


3.256427(2) 
-1.398(15) 



Table 1. Total magnetic moments (in (in) of 7 Be and 7 Li for 
different choices of relative angular momenta compared with 
experimental value (last column). Even forgetting parity con- 
siderations that would of course leave only the I = 1 state 
as a possible candidate, it is clear that the best agreement is 
obtained in both cases by considering a relative p-wave. 



We calculate the relative motion wavefunction of the 
ground state, first excited state and continuum states by 
solving numerically the radial Schrodinger equation with 
the intercluster potential given in (jl . 1|) . In order to yield 
the resonant 7/2~ and 5/2 _ states, we have to use a dif- 
ferent set of parameters for the Woods-Saxon (Vws = 
-64.68 MeV) and spin-orbit (V; s = 1.94 MeV) poten- 
tials, but keeping the same geometry. With this choice 
we find that the width of the and f$/2 resonances 
are in reasonable agreement with the experimental data 
(see table ©). A number of ground state properties are 
found to be in good agreement with measured observ- 
ables. We also calculate the B(E2,3/2~ -> l/2~) and 
B(Ml,3/2~ — > l/2~) transition strengths between the 
ground state and the first excited bound state. All the cal- 
culated observables, obtained using the expressions given 
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Quantity 


this work 


experiments 


(r z ) A+B (fm) 


2.48 


2.48 ± 0.03 


{tYa+b (fm) 


2.52 


2.52 ± 0.03 


gmatter ( fm ^ 


-11.83 




QCharge 


- 4.79 




A* (Mat) 


-1.53 


-1.398 ± 0.015 


B(K2,3/2~ -> 1/2") (e 2 /m 4 ) 


18.3 




S(Ml,3/2- ->■ 1/2") (a* 2 ) 


1.86 


1.87 ± 0.25 


r(7/2~) (keV) 


~ 90 


175 ± 7 


r(5/2 - ) (MeV) 


~ 0.8 


1.2 



Table 2. Static properties of 7 Be in a dicluster model and elec- 
tromagnetic transition strengths between bound states. Defini- 
tions of some of these observables can be found in the appendix. 
Experimental values are taken from |24l ; [2BI ; I27I ; I23T ) 




E (McV) 



in appendix, are collected in table ([2]) with a comparison 
with known experimental data. 

As it has been done in Ref. (@) for 7 Li, besides the qua- 
drupole transition to the first excited bound state given in 
Table (J5J), we also investigate the electromagnetic response 
of 7 Be to continuum states. Figure ([l} shows differential 
reduced transition probability (see Ref.|28|)), dB(El)/dE c , 
for continuum states allowed by an El transition. In this 
case there are no resonances in the low-lying continuum 
and the visible peaks have a non-resonant nature (fl3l). 
Dotted line refers to continuum s x y 2 state, dashed line 
and long-dashed refers to d 3 / 2 an d d 5 / 2 states respectively. 
The solid line gives the sum of all contributions. The total 
integrated non-energy weighted B(E1) is 0.226 e 2 fm 2 (See 
next section for comments on sum rules). 

Figure shows differential dB(E2)/dE c for contin- 
uum states allowed by an E2 transition. Long-dashed and 
dot-dashed lines refer to f$/2 and fi/ 2 states respectively, 
while dotted and dashed lines refer to P1/2 and p 3 / 2 states 
respectively. The / states take contributions both from the 
resonances and from the non-resonant part of the contin- 
uum, while the peaks of the p states arise purely from non 
resonant transitions. We show also in the inset the full f-j / 2 
resonance on a different scale because it reaches a maxi- 
mum of about 150 e 2 fm 4 /MeV. Finally solid line gives the 
sum of all contributions. The total integrated non-energy 
weighted B(E2) is 113.5 e 2 fm 4 to be compared with the 
B(E2,3/2~ -> l/2~) value of 18.3 e 2 fm 4 . Non resonant 
continuum states P1/2 and P3/ 2 are dominant at very low 
continuum energy, while they become negligible at higher 
energies around and beyond the f§/ 2 peak. 

In addition to electric excitations we also consider mag- 
netic ones. Magnetic dipole interactions give normally con- 
tributions to the transition probabilities that are almost 
comparable with electric quadrupole and smaller than elec- 
tric dipole. It appears however, already in old calcula- 
tions (fl4r ) for the photo-disintegration of the deuteron (the 
smallest dicluster nucleus), that they are relevant at as- 
trophysical energies (say below 200 keV) where they even 
turn out to dominate over the electric dipole. For this 
reason we decided to investigate this theme to better un- 



Figure 1. Differential dB(El)/dE c in e 2 fm 2 /MeV for transi- 
tions from the ground-state to the continuum. Energies are in 
MeV, referred to the threshold for breaking into the a— 3 He 
channel. Separated contributions are indicated in the legend. 
Solid line shows the sum of all the electric dipole contributions. 
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Figure 2. Differential dB(E2)/E c in e 2 fm 4 /MeV for transi- 
tions from the ground-state to the continuum. Energies are in 
MeV, referred to the threshold for breaking into the a— 3 He 
channel. Separated contributions are indicated in the legend. 
Solid line shows the sum of all the electric quadrupole contri- 
butions. 



derstand the role they might play in low-energy nuclear 
reactions involving dicluster systems (cfr Sect. 4). 

The expression for the magnetic operator can be ob- 
tained from the multipoles expansions of the electromag- 
netic field with sources as (cfr. Ref. (fl5l). ch.17): 



M(Ml)=J2lV(r?Y x (d,<j>))} 



<ffi h 



-rriiec(\ + 1) 



where the sum is over all particles (with cordinates r, 9, <f>), 
each with an effective charge ^ff-i ma g ne tic moment fii 
and mass m^. As usual A indicate the multipolarity, e is 
the electron charge. 
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In the present case the wave functions are not ex- 
pressed in a single-particle proton or neutron basis, but 
rather in a two-cluster basis with the a (which is spinless) 
and either the triton or the helium-3. Therefore we can- 
not adopt the same simplifications as in Ref. dH), but we 
have to rewrite the magnetic operator in a cluster form. 
For the dipole case this is not difficult and the magnetic 
operator can be recast in the form: 



M(M1) = VSY x + (j, 2 + UnGL) 
= V3Y x ( m GJ + (2/i d - G)S) 



(1.3) 



where the summation over the A particles has been sim- 
plified by partitioning the protons and neutrons into two 
spherical clusters (see appendix). Here Y indicates vec- 
tor spherical harmonics, /zi and /i2 are the intrinsic static 
magnetic operator of the clusters (in nuclear magnetons, 
Hn) and L = J — S is the cluster-cluster orbital angular 
momentum operator. The total angular momentum and 
spin of the cluster are indicated with J and S respectively. 

Z -\- Z 

The coefficient G = amaI acts 

as an effective cluster- 
cluster orbital gyromagnetic factor. Since the a has null 
spin and null magnetic moment, only the intrinsic mag- 
netic moment of the other cluster appears (ji c i = ^2)- 

By calculating the reduced matrix elements along the 
lines of Ref. (|15l). one gets: 

B{Ml) = 9(2jV + 1) | (lfj f | k,3i) I 2 (if\\Y m 2 - 



fiN G(-l) t '+'' + i'- 1 {i i ji\\J\\ei,j 




+ (2fi cl ~ f j, N G)(s f \\s\\ 



where li,ji,lf and jf are initial and final orbital and total 
angular momenta. The selection rule given by the overlap 
and the reduced matrix element of the vector spherical 
harmonic imply that one can only populate states with 
£i — if and ji ^ jf. In our specific case, where the initial 
state has tj — P3/2, we can only reach states with P1/2 
character. The corresponding B(M1) is : 



B(M1) 



4tt 



(2n cl - m G) 2 I (l f jf I kJi) | 2 (1 . 5) 



that coincides to the formula given in Refs. (dH). The cal- 
culation of the B(M1; 3/2 — > 1/2) to the bound state (that 
is within the two isospin partner levels) gives an excel- 
lent agreement with experimental measurements as seen 
from table [H Figure §5§ shows continuum states excited 
in Ml transitions for 7 Be. Only transitions to the pi/2 
states are allowed because of orthogonality between the 
ground-state P3/2 and continuum p 3 / 2 (see the next sec- 
tion for a more detailed discussion). The total non-energy 
weighted B(M1) to the continuum states, integrated up 
to 50 MeV, is about 0.00408 [i 2 , much smaller than the 
magnetic dipole strength to the first bound state. 



B(MI)=1.86 



p3/2->pl/2 



4 5 
E (MeV) 



Figure 3. Differential dB(Ml)/dE c in /x 2 /MeV for transitions 
from ground-state to the continuum. Energies are in MeV, re- 
ferred to the threshold for breaking into the a— 3 He channel. 



1.2 7 Li 

In a the previous paper we have investigated dipole 
and quadrupole electric transitions of the nucleus 7 Li. We 
refer here to the same parameter sets for the potential and 
we improve our calculations by adding results for magnetic 
dipole transitions to the continuum for the sake of com- 
pleteness. We note that 7 Li and 7 Be have the same qual- 
itative behaviour because of their similar internal struc- 
ture. An estimate of static magnetic dipole moment of 
7 Li ground-state obtained by using formula (|5 . 20|) gives 
: nij Li) — 3.37/iAr that is in good agreement with the ex- 
perimental value 3.256^jv given in (|26[ ). Again this fact 
confirms that in a simplified dicluster picture the I = 1 
relative angular momentum gives a good magnetic dipole 
moment as in 7 Be. 

7 Li and 7 Be have also the same qualitative behaviour 
with respect to transitions to the continuum: compare 
plots in Ref. (Q) with figures ([T]) and {2} in the present 
paper and compare figure |3]) with figure ((4]) . In particular 
the last one shows the magnetic dipole transition strength 
for transitions between the ground-state and the contin- 
uum. Only transition to P1/2 continuum states are allowed 
because of orthogonality between the ground-state p 3 / 2 
and the P3/2 continuum. This point was discussed, for an- 
other system, also in Ref. (flU). 

The total integrated B(M1) strength amounts to 0.0029^i 2 
again much smaller than the B(M1) to the first excited 
bound state. The magnetic dipole contribution is, in both 
nuclides, not as important as for the deuteron. 



2 Molecular sum rules 

Useful tools in nuclear physics are sum rules (fl7l : list [H; 
I20I ). In particular energy weighted sum rules (EWSM) for 
an electromagnetic interaction can be expressed in terms 
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p3/2->pl/2 




3 4 5 
E (McV) 



Figure 4. Differential dB(Ml) /dE c in fi 2 /MeV for transitions 
from ground-state to continuum. Energies are in MeV, referred 
to the threshold for breaking into the a — t channel. 

of proper expectation values in the ground state 



S(ttA) = ^(Ef - Ei)B(nX; i -> /) 
/ 

= 1(0 | [6(wX),[H,6(ttX)}] |0) 



(2.6) 



where i stands for the initial state and the sum is extended 
over all reachable final states /. Here H is the hamiltonian 
operator of the system. The energies of the initial and 
final states are Ei and Ef and the operator 0(irX) in the 
double commutator is the one that induces the transition. 
For electric dipole and quadrupole interactions, equation 
(|2 . 6p leads respectively to (|17l ): 



S(El) = 
S(E2) = 



9 h 2 NZ 



Air 2m 
47T 2m 



A 

Ze 2 {r 2 ) ch 



(2.7) 



where m is the nucleon mass and (r 2 ) ch is the mean square 



charge radius. Our calculations show that the EWSR (|2 . 7p 
are clearly not exhausted by considering only the transi- 
tions which imply changes in the relative motion (see table 
HP). We assume that the suitable sum rules in this case 
are the energy weighted molecular sum rules (EWMSR) 
that can be obtained by removing the contributions of the 
individual clusters (|18f h If a nucleus with mass and charge 
(A, Z) is split into two clusters (Ai, Zi) and (A2, Z2) then 
the EWMSR is in general defined as 

S mo i(EX,A 1 +A 2 ) = 

= S{EX, A x + At) - S(EX, Ax) - S(EX, A 2 ) . (2.8) 

For electric dipole (flU ) and quadrupole interactions one 
obtains: 

Smol(El) 

S mo i(E2) 



9 (Zt_A 2 




The EWMSR for the quadrupole case is obtained by 
using formula (|5 . 18[) into eq. (|2 . 8|) . Note that in Ref. 
(|18r ) the expectation value of (R) 2 was approximated, in a 
phenomenological fashion, by the intercluster equilibrium 
distance, So- 

Energy weighted sum rules for the magnetic dipole (f2lh 
and higher multipolarities (|22h have been evaluated in a 
pure single-particle scheme with a shell model hamiltonian 
containing a spin-orbit part. The magnitude of the spin- 
orbit splitting is directly proportional to the total energy- 
weighted strength. In our cluster picture, one can follow 
a similar idea to evaluate the Energy- Weighted Molecular 
Sum Rule for the magnetic dipole interaction. The model 
hamiltonian is H = Hq + V so (R)L ■ S, where Ho contains 
the kinetic term and the central part of the potential, 
while the operator that promotes magnetic dipole excita- 
tions amounts to the third component of O = J2 i=1 2 (GLi-{ 
2/XiSj), where G (defined above and also in 15 . 20[) con- 
tains the combined effects of effective charges and reduced 
mass. The expectation value in the ground state of double 
commutator gives 



SmoliM^Ax 



A 2 ) = -{0\[O z ,[H,O z }}\0) 



= --(G-2 f i 2 ) 2 {0\V so (R)L-S\0) 



(2.9) 



We have calculated all the contributions of El and 
E2 transitions, including also the contribution of "vir- 
tual" states, or states that are not allowed by the Wilder- 
muth's rule, (Osi/2, lsi/2> 0^5/2i 0d 3 /2 for electric dipole 
and 0p3/2,0pi/2 for quadrupole). For magnetic dipole in- 
teraction we have to consider the virtual state 0pi/ 2 , the 
first excited state lpi/2 and the continuum P1/2, but we 
do not need to consider the virtual state 0p3/ 2 . The rea- 
son is the following: in the reduced matrix element for an 
MA interaction the radial contribution is (f\r x ~ 1 \i). For 
magnetic dipole interaction this contribution is simply the 
overlap between the initial state, i, and the final state, /. 
If the states carry the same angular momentum and spin 
quantum numbers then their overlap is zero by orthogo- 
nality, because they are different eigenstates of the same 
hamiltonian (spin-orbit potential included). For the same 
reason the virtual state 0pi/ 2 does indeed contribute: since 
the spin-orbit part of the potential is different, the two 
hamiltonians are different and one cannot use the prop- 
erty of orthogonality, but rather has to calculate the actual 
overlap. 

In table © we compare EWSR and EWMSR with the 
values obtained in our model. We find that the low-lying 
dipole strength exhausts almost entirely the EWMSM, but 
represents only a small fraction of the EWSR. In the same 
way the quadrupole strength coming from changes in the 
relative motion of dicluster configurations exhausts a large 
fraction of the EWMSR, but only a small part of the stan- 
dard EWSR. 

The energy weighted molecular sum rule for magnetic 
dipole transitions is fulfilled quite well by our calculations: 
% in the case of lithium and 94% in the case of berillium. 
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'Li 


'Be 


EM A 


EWMSR 


EWSR 


ours 


EWMSR 


EWSR 


ours 


El 


1.02 


36.7 


1.01 


1.02 


36.7 


1.00 


E2 


1120 


2105 


481.6 


1424 


3018 


639 


Ml 


1.36 




1.20 


0.87 




0.82 



Table 3. Comparison between EWSR, EWMSR and our calculations for 7 Li and 7 Be. Values are in MeV 2 fm 3 for El, MeV 2 fm 5 
for E2 and fj, 2 N MeV for Ml. 



Practically all the contribution comes from the excitation 
to the first excited bound state. 



3 Photo-dissociation and radiative capture 

The knowledge of the electromagnetic response to contin- 
uum is a basic ingredient to describe break-up processes. 
In fact, in kinematic conditions where the process is domi- 
nated by the Coulomb interaction (for example at very for- 
ward angles), the break-up probabilities become directly 
proportional to the B(EX) values. Independently from 
kinematic conditions this is also the case for two other 
processes of fundamental importance for astrophysics that 
can be therefore studied within the present model: the 
photo-dissociation and its inverse process, the radiative 
capture. 

The electric and magnetic response functions obtained 
so far can be used to calculate the cross sections for these 
two processes. The radiative capture cross-section for type 
(n = Electric or Magnetic) and multipolarity A, can be 
expressed (|28[ ) as 



Ccapi(7rA, Ec) 



2(2j + 1) 



Vph.dis(K\, E c ) 



(2j a + l)(2jci + 1) \k a - cl 

(3 . 10) 

where fc a _ c i is the wave-number for the relative motion of 
the two clusters (one is always an a particle in the present 



case) and 



2(2j + l) 



is a spin factor from the detailed 



(2i„+i)(2i ci +r) 

balance principle. This expression relies on the knowledge 
of the photo-dissociation cross-section that may be ex- 
pressed as follows: 



Tph.dis{^^,E C ) = (27r)"' 



(A + l) 



k 

[A(2A + 1)!!] 2 7 



2A- 



dE c 
(3.H) 

where dB is the reduced transition probability (de- 

fined as in Ref. |28l)). fc 7 is the photon wave-number, A is 
the multipolarity of the transition, E c = h 2 k 2 t _ cl /2^ r is 
the relative energy of the two clusters in the continuum 
and E-y = E c + Eb is the photon energy. Here we have 
used the notation /i r for the reduced mass and Eb for the 
binding energy. 

To evaluate the total photo-dissociation cross-section, 
starting from a given bound state, one has to sum over all 
possible angular momentum transfer and possible allowed 



spins of the final states in the continuum. In the calcu- 
lation of the total capture cross-section, in addition, one 
has to include all the bound states. 

Fig. ([5]) shows a schematic picture for the two processes 
(left side), and the respective cross sections (right side). 
They are very similar at a qualitative level, but the differ- 
ent binding energies of the ground states (mostly due to 
the different charge) makes up for changes in the position 
of the resonances and in the intensity of the transition to 
the low-lyng continuum. 

This is evident from the comparison of fig. [T] of the 
present paper and fig. 1 of our older work on 7 Li (|9|). The 
BE(1) to low-lying continuum states has slightly different 
contribution from the s and d states that explains the 
differences in the cross-section at small energies. Indeed it 
is known that the shape and position of the maxima in the 
response are affected by binding energy, initial and final 
angular momenta and Coulomb barrier height (|28l: l29h. 

Electric dipole clearly dominates the cross-section, al- 
though the electric quadrupole shows peaks due to the 
presence of two low-lying resonances. The electric qua- 
drupole and the magnetic dipole are of comparable size 
and are nevertheless not negligibly small in comparison 
with the other contributions. 

The photo-dissociation cross-section is shown only from 
the ground state, but the radiative capture can go directly 
to the ground state or proceed through the first excited 
bound state. In the latter case, due to the small energy 
difference, the cross sections to the P1/2 and p 3 / 2 are prac- 
tically proportional and the proportionality is fixed by the 
phase space factor depending on the total spin, j. 



4 S factors 

In the stellar environment the occurrence of high temper- 
atures and the large density of light isotopes makes the 
a— capture a very likely process. Tritium and 3 He partici- 
pate in (a, 7) reactions forming mass A = 7 isobars, which 
are one possible getaway toward heavier nuclei. It is there- 
fore important to understand well these processes and to 
have simple models to treat the reactions. Once again the 
dicluster picture qualitatively and quantitatively entails 
the correct physics ingredients: indeed the astrophysical 
S- factor (a reparameterization of the cross-section, defined 
by S = e re i<T cap t(e re i)e 27rri , where rj is the Sommerfeld pa- 
rameter) for the reaction 3 He(a,7) 7 Be is reasonably well 
reproduced as can be seen in figures ([fJ]). Good results 
are also obtained for the partner reaction, 3 H(o;,7) 7 Li, 
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' 7 3 

Li + y -> a + t Be +y -> a + He 




E c (MeV) E c (MeV) 

Figure 5. Photo-dissociation and radiative capture processes (left pictures) and cross-sections (right) plotted against the 
continuum energy for reactions involving 7 Li and 7 Be. The E1,E2 and Ml contributions are shown separately and the total 
practically coincides with the electric dipole contribution. The E2 resonances are seen in each picture. The photo-dissociation 
is shown only from the ground state, but the radiative capture can go directly to the ground state or proceed through the first 
excited state. 
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Figure 8. Integrand of the response function for the transition 
S1/2 — * P3/2 in 7 Be for a given energy E c = 50 keV (upper 
panel) and S-factor for the same transition in the low-energy 
regime. The two curves are calculated including (dashed) or 
excluding (solid) the nuclear potential in the evaluation of the 
continuum wavefunction. 



as shown in 0. Each curve has both a positive and a 
negative aspect: while the result for Be has a good magni- 
tude, but a non-optimal slope, the result for Li has a good 
shape, but it overestimates the data. Our results do not 
contain any rescaling of the S-factor or, in other words, 
we assume a spectroscopic factor for the a+cluster con- 
figuration equal to one. This might be the reason for the 
observed deviation in the case of 7 Li, where the 6 Li+n 
component may affect the total magnitude although pos- 
sibly not the overall shape. The reason for the flat shape 
of the 7 Be curve might also reside in the prescription we 
have used to treat the dicluster system. Perhaps it is an 
indication that other components are affecting the prop- 
erties at low energies although the results in table II were 
very satisfactory. 

The various sets of experimental measurements that 
are reported in these figures together with the calculated 
values, are taken from the NACRE database (|30t ). where 
an extensive list of references can be found. A simple 
parabolic fit of the model curves give a S(0) of 0.42 and 
0.13 keV b for beryllium and lithium respectively. 

Other more refined approaches, like the R- matrix cal- 
culations (see discussion in Ref. (112) )■ or the Multichan- 
nel Algebraic Scattering approach (fill ) have been applied 
with success to the same reaction. Although our model is 
much simpler than the cited approaches, it has a compa- 
rable predictive power. Of course this might still not be 
enough for elaborated nucleosynthesis models that must 
consider a network of nuclear reactions with very precise 
fits. In that case a more precise determination of the S- 
factor might be necessary and one must go beyond the 
simple cluster picture. 

The whole idea of parametrizing the capture cross- 
section in terms of the S-factors is precisely to smooth 
out the effects of the Coulomb interaction (by means of 
the exponential term containing the Sommerfeld param- 



eter). The important effect of the nuclear interaction in 
the low-energy regime can be seen in Fig. |8] and cannot 
be neglected. The upper panel displays the integrand of 
the response function as a function of r for a given en- 
ergy of E c = 50 keV in the continuum. We have selected 
the S]/2 ~~ * P3/2 transition in 7 Be for the sake of simplic- 
ity. Clearly the nuclear potential does not affect the tail 
of the wavefunction, that matches a pure Coulomb wave- 
function at large distances, but at distances smaller than 
the position of the Coulomb barrier, the presence of the 
nuclear potential affect the continuum wavefunctions and 
therefore the response in the continuum. The S-factor, in 
the low-energy regime, is displayed in the lower panel. 
The nuclear potential is seen to lower the S-factor sig- 
nificantly and cannot be excluded from calculations, as 
already found by several authors. 



5 Conclusions 

An elementary dicluster model is found to provide an ef- 
fective description of static as well as dynamic properties 
of some light nuclei and can be used with confidence in re- 
action models. Although ab initio many-body models are 
expected to provide better results in the general case, we 
have shown that, for example, the astrophysical S-factor, 
which is a crucial quantity for interdisciplinary applica- 
tions of nuclear physics, is reasonably well reproduced by 
the dicluster model. 

In particular we have investigated the magnetic prop- 
erties of A = 7 isobars, introducing a new molecular mag- 
netic dipole sum rule, that has a close analogy with the 
single particle magnetic sum rule (j22l ). and finding that 
(as expected) the contribution to photo-dissociation and 
radiative capture processes is rather small. This is in con- 
trast with the strong Ml photo-dissociation peak at low 
energies in the case of deuteron, that has motivated our 
investigation. 



Appendix: derivation of standard formulae for 
dicluster nuclei 

We summarize here a few important formulas valid for 
dicluster nuclei (fH: I27I). They are often used in the litera- 
ture, but the derivation, albeit elementary in some cases, 
is rarely reported. 

Matter radius 

The mean square matter radius of a nucleus with mass 
number A is defined as 



A 



(5 . 12) 



Assuming that the nucleus made up by two clusters with 
mass number Ai, A 2 and A=Ai + A 2 , Eq. (|5 . 12[) be- 
comes 



" A = A 



(5 . 13) 



A.Mason, et al: Response to the continuum for A — 7 isobars 



9 



0.8 



0.6 



J3 
> 



U 0.4 



0.2 






calculated 




o 


PA63 




□ 


H059 




o 


HI88 




V 


AL84 






KR82 






NA69 




< 


OS82 




> 


R083 




A 


LU07 






3 7 

He (a,y) Be 



_L 



_L 



0.5 



1.5 



E c (MeV) 

Figure 6. S-factor for the 3 He(a,7) 7 Be reaction. 




calculated 

GR61 

SC67 

BU87 

BR94 



Figure 7. S-factor for the t(a,^) 7 \A reaction. 



where Rij are position vectors measured from the center By inserting (|5 . 14|) in (|5 . 13|) and upon taking the ex- 
of mass of each cluster, A1.2 respectively; R is the inter- pectation values rather than the operators, one gets 
cluster distance that can be split into two vectors, Ra ± 
and Ra 2 related by 



Ra ± 
Ra 2 



A 



R 



Ai 
A 



R 



(5 . 14) ^) Al+A2 = ^) Al + ^) A2 + ^^) (5 . 15) 
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Charge radius 

The mean square charge radius of a nucleus with atomic 
number Z is 

Ez 2 



' ch 



(5 . 16) 



For a nucleus composed by two clusters with atomic num- 
bers Zi, Z2 and Z=Zi+Z2 eq. (|5 . 16[) becomes 



' ch 



(5 . 17) 



Using (|5 . 14[) and taking the expectation values of oper- 
ators, one finds 



„2\ch 



J l i2\ch 



J2 lJ2\ch 



/A1+A2 



z 



(R 2 



Z 



(5 . 18) 



Magnetic dipole moment 



The magnetic dipole moment operator is defined as 



A+B 
/2 W 



A+B 



(5 . 19) 



i=l 



with g l = 1 for protons, g l = for neutrons, g s = 5.58 for 
protons and 17 s = —3.83 for neutrons. The orbital angular 
momenta are defined as U = riXpi, where r< = Ra^ +R-i 
for cluster Ai and Tj = Ra 2 +Rj i0r cluster A 2 . The same 
rule can be applied to momenta pi. Labelling with /xi, fi 2 
the intrinsic magnetic dipole moments of each of the two 
clusters, (|5 . 19|) becomes 



/x = Hi + fi 2 + 



ZiAl 



Z 2 A\ 



AA 1 A 2 



(5 . 20) 



where one has isolated in the third term the contribution 
from the relative orbital angular momentum, L. The to- 
tal magnetic moment is therefore simply the sum of the 
magnetic moments of the two subsystems plus the con- 
tribution from the cluster-cluster orbital motion, which 
vanishes if the two cluster are in a relative s-state. 



and, by using the definition of the spherical harmonic and 
by splitting the sum over the total nucleon number into 
two sum over the nucleons of each cluster, we obtain two 
relationships that contain the intrinsic quadrupole mo- 
ments of the clusters and a term that expresses the qua- 
drupole operator for the relative motion, namely: 



qmat = gmat + , 



A1A 2 . + A 2 A\ 
A 2 

12 , rj_ A2 



R 2 Y 20 (©,#) 



Q ch = Qt + Qt + ZiA2 + 2 Z2Ai R 2 y 20 (o,<p) 

(5 . 23) 



With these operators one can calculate the quadrupole 
moments by evaluating the expectation value in the ground 
state with the maximally aligned magnetic substate: 



16tt 



(L, M = L | Q | L, M = L) (5 . 24) 



For example in the case of the matter quadrupole moment 
one has: 



r\mat.mo rymat.mo , ryraat.rao 



A\A\ 



A 2 A\ 



A 2 



(i? 2 )2(2L + l)(-l) i 



L 2L 
-L0L 



L2L 


(5 . 25) 

where L is the angular momentum of the relative motion. 
In the case of the two mirror isobars lithium and beryl- 
lium, we have L = 1 and therefore the above formula 
reduces to 



Q' 



ryn\at.n\o 



r^mat.mo 
^A 2 



2 A X A\ + A 2 A 



i(i? 2 



5 A 2 

(5 . 26) 

Analogous expressions may be derived for the charge qua- 
drupole moments that differ only in the coefficient that 
depends linearly on charges rather than masses. Notice 
also that, while the coefficient of the relative motion part 
for the charge quadrupole moment is different for lithium 
and beryllium, the one for the matter quadrupole moment 
is the same. Another difference comes from the average 
square radius of the relative motion wave function, which 
is not identical for the two nuclei. 
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